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Abstract 

The autocorrelation function, A{t)^ measures the overlap (in Hilbert space) of a time-dependent 
quantum mechanical wave function, 'ip{x,t), with its initial value, ^{x,0). It finds extensive use in 
the theoretical analysis and experimental measurement of such phenomena as quantum wave packet 
revivals. We evaluate explicit expressions for the autocorrelation function for time-dependent 
Gaussian solutions of the Schrodinger equation corresponding to the cases of a free particle, a 
particle undergoing uniform acceleration, a particle in a harmonic oscillator potential, and a system 
corresponding to an unstable equilibrium (the so-called 'inverted' oscillator.) We emphasize the 
importance of momentum-space methods where such calculations are often more straightforwardly 
realized, as well as stressing their role in providing complementary information to results obtained 
using position-space wavefunctions. 
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I. INTRODUCTION 



Despite widely varying approaches to the subject, every introductory quantum mechan- 
ics text of which we are aware discusses the subject of wave function normahzation and 
conservation of probabihty, proving that 

{i^Mt) = (^ol^o) = 1 (1) 

provided the Hamihonian describing the system is Hermitian. This important relationship 
is perhaps most often discussed in position-space using 



{^t\^t) = I r{x,t)^{x,t)dx = l (2) 
but, using the Fourier transform. 



hoc 



V2TihJ- 

it is equally well written in the form 



oo 



{^t\^t)= / <P*M<P{p,t)dp=l. (4) 

J — oo 

For bound state problems, where any time-dependent state can be expanded in terms of 
energy eigenstates, with quantized energies, we can also write 

/ + 00 
[un{x)]* ^{x,Q)dx (5) 

SO that the normalization condition is given by 

oo 



n=l 



An important related concept, namely the autocorrelation function A{t), is defined 
as the overlap of a time-dependent state, ifji, with its initial value, ■i/'o, namely 

/■+00 r+oo 

A(t) = (^il^o) = / r{x,t)^{x,^)dx= / 0*(p,t)</>(p,O)cip (7) 



and is extensively used in the study of the time-development of quantum wave packets in 
bound state systems. One clearly has A(0) = 1, by definition, if the initial wave function is 
itself properly normalized, but generally \A{t)\ < 1 for later times as the wave packet de- 
velops in time and the different energy /momentum components evolve differently. Besides 



being of obvious theoretical value in the analysis of time-dependent systems, the autocorre- 
lation function is physically important because it is very directly related to the observable 
ionization signal in the pump-probe type experiments where the time-development of atomic 
wave packets is studied experimentally 

a. a. 

For a single (bound) energy eigenstate or stationary state, where ^{x,t) = Un{x)e~^^"*/^^ 
one clearly has 

A{t) = e+*^"*/^ giving \A{t)\^ = 1 . (8) 

For more general bound state systems as in Eqn. (0), however, the most appropriate form 
is 

oo 

A(t) = (^,|^o) = 5^|anre+^^"*/'^ (9) 

n=l 

SO that information on the non-trivial time-development of the wave packet is encoded in the 
energy eigenvalue spectrum. For wave packets which are constructed from energy eigenstates 
centered around some large value of n = uq, one can write 

E{n) ^ E{no) + E'{no){n - Uq) + ^-^{n - no)' + • • • (10) 

which gives the time-dependence of each individual quantum eigenstate as 

= exp (— if2o^ ~ 27ri(n — nQ)t/Tci — 2ni(n — noft/Trev + ■ ■ ■ ) (H) 



E(no)t + (n - no)E'(no)t + ^(n - nQfE'\no)t + 



in terms of which the classical period and quantum mechanical revival time (discussed below) 
are given respectively by 

Tri = — ; — TT and Irev = ■; ; — m- (12) 

\E'{no)\ \E"ino)\/2 ^ ^ 

and the common first term, exp(— z^ot) = exp{—iE{nQ)t/h), is an unimportant overall 
phase. The most familiar example is the harmonic oscillator, with En = {n + l/2)huj, and 
any wavepacket in this system is periodic with period Td = 2n/u and Trev — > oo since 
E"{no) = 0. 

This form of the autocorrelation function is especially useful in the context of quantum 
wave packet revivals l|, that is, systems where initially localized states which have a short- 
term, quasi-classical time evolution, can spread significantly over several orbits, only to 
reform later in the form of a quantum revival in which the spreading reverses itself, the 



wave packet relocalizes, and the semi-classical periodicity is once again evident. The presence 
of an approximate quantum revival at a later time, t = Trev, is indicated by |74(Tre«)| ~ 1, 
accompanied by evidence of the return of the short-term periodicity in A{t) with the classical 
period, T^. The phase-structure of A{t) can also yield useful information, as parametric 
plots of -Re[A(i(:)] versus Im[A{t)] (Argand diagrams) can provide striking visualizations [3] 
of the highly correlated Schrodinger cat-like states which evolve at fractional multiples of 
the revival time, the mathematics of which was first worked out in detail by Averbukh and 
Perelman P]. 

Such revival phenomena have been observed in a wide variety of physical systems, es- 



pecially in Rydberg atoms p - j8[. In many cases, the autocorrelation function is closely 
related to experimentally observable quantities , , and its use has become very familiar 
for analyzing model systems exhibiting exact or approximate wave packet revivals. 

We note that early investigators commented on the existence of a 'quantum recurrence 
theorem' which utilized the notion of 'distance' (in Hilbert space) between a time- 
dependent quantum state and its initial value as 

r+oo 



oo 



2 1-3? 



x,t)\'^dx+ / \^{x,0)\'^dx 

J — oo 

oo f+OO 

tp*{x,t)tp{x,0)dx- / ij*{x,0)tp{x,t)dx (13) 

J — oo 

tlj*{x, t) tpix, 0) dx] ) = 2(1 - ^[A{t)]) 



which is also related to A(t); we also note that Baltz has considered similar ideas in a 
more pedagogical context. 

In contrast to the research literature, where A{t) is now a standard tool, there are few, 
if any, examples of the evaluation of the autocorrelation function for the many familiar 
and frequently studied model systems of introductory quantum mechanics. Besides being 
useful as a diagnostic for the rate of time-evolution of a quantum state, such calculations of 
A{t) can also help answer questions such as How similar are two quantum states?, both in 
magnitude and phase, and we will stress the complementary roles that position-space and 
momentum-space approaches can have in addressing such issues in what follows. 

In this note, we will use the special properties of Gaussian wave packet solutions to 
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explicitly evaluate A{t), obtaining closed form expressions, for four familiar and accessible 
model systems. We focus on the cases of a free particle (Sec. II), a particle undergoing 
uniform acceleration (Sec. Ill), a particle in a harmonic oscillator potential (Sec. IV), and a 
system corresponding to an unstable equilibrium (the 'inverted' oscillator, in Sec. V.) The 
first two cases do not correspond to bound state systems, but do provide useful results for 
comparison to the short-term time- dependence of wave packets in systems such as the infinite 
well 0], 11| - 17 1 and the so-called 'quantum bouncer' - 21 1 where, for small times at 



least, the wave packet propagation is similar to the corresponding unbound case. For the case 
of the harmonic oscillator, any wave packet solution (Gaussian or not) is explicitly periodic 
with Tci = 27t/uj and the wavepacket never enters a truly 'collapsed' phase; the oscillator 
does, however, provide a useful explicit example illustrating the exactly periodic behavior 
of A{t), seen more approximately in many bound state systems. For example, Nauenberg 
^ has provided an elegant general description of the medium-term time-development of 
A(t) for a general, one-dimensional, bound state system and the exact results presented 
here for the harmonic oscillator can be used as an efficient 'benchmark' for comparison to 
that more general analysis. While we will focus on closed-form results for oscillator wave 
packets (obtained using propagator methods), we can also make contact with the expansion 
in eigenstates in Eqn. as well. 

Such studies of the general behavior of A{t) for many standard ex amp les are also useful as 



: amp ies are 

H 



they con,ple„.ent exiting woA on the .ate ot wave packet .pteading ^21 - M and e.pectal,y 

on the time evolution of quantum states |2q], [2a|. In this last context, the example provided 
here for the free particle (in Sec. II) can be used as a specific case to confirm a (hopefully) 
well-known result for isolated quantum systems js^] - j^]; namely 

K^,|^o)r>cos^(^) for 0<t<^ (14) 

where AH = a/ {H'^) — {H)"^ is the uncertainty in the free-particle energy of the wave packet. 

Because our presentation here will use the same notation and many of the same methods 
as the companion paper js^, we will refer extensively to results from that paper, especially 
for properties of the standard Gaussian wave packets we utilize. 
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II. FREE-PARTICLE GAUSSIAN WAVE PACKETS 



The most general free-particle, momentum-space and position-space Gaussian wave pack- 
ets, with arbitrary initial values of (a;)o = Xq and (p)o = Po, can be written [s^ in the form 



-a2{p-po)^/2 -ipxo/h -ip'^t/2mh 



0) 



a 



g-a2(p_pQ)2/2 ^-ipxo/h 



and 



a/ y/7Tah{l + it/to 
1 



^ipo(x-xo)/h -ip'^t/2mh -(x-XQ-pot/mf' /2{aKf(l+it/t(,) 



Tiah 



^ipo{x-xo)/h ^-{x-xq)'^ /2(ah)'^ 



(15) 
(16) 

(17) 
(18) 



where to = TTiha^ defines the spreading time. The calculation of A{t) is done most straight- 
forwardly in momentum-space where 



A{t) 



a 



exp 



ia^plt 



and the modulus-squared is then given by 
\A{t)f ^ 



_2to{l-it/2to) 



Vl + {t/2to) 



exp 



2 2 (^2^0 

-2a Po 



l + (t/2to)^)J 



The same result can, of course, be obtained in position-space, with 



A{t) 



ip*{x, t) ip{x, 0) dx 



^ip^t/2mh 



exp 



ipot/m) 



4:{ah)^{l-it/2to) 



(19) 



(20) 



(21) 



which is easily seen to be identical to Eqn. (jl9p with a minimum of manipulation. 

We note that for times satisfying t << to there is an increasing exponential suppression 
of the overlap between ipt and ipo, but the exponential factor does 'saturate' for long times, 
giving 

9±. r r,2 1 1 

since Apo = — ^ . (22) 



\Ait»2to)\'—*-^exp 



a 



V2 



The asymptotic form of the exponential factor can be understood by noting that the 'distance 
in position space' between the initial 'peak' at (x)o = xq, and that at later times when 
{x)t = Xq +pot/m, grows linearly with t, while for long times the position spread. 



Axt = AxoVl + (t/to)^ 
6 



^ t 
Axo— , 
to 



(23) 



grows in the same way. This leads to factors in the exponent of the form 

{x{t)-x{0)f {pot/mf [ poto ^2 . A /^.x 

~ ' ' ~ (Pott) smce Axo = ■ (24) 



This argument is most transparent using the position-space wave functions where the 
exponentially small overlap of magnitudes is clear, but the same suppression arises in the 
momentum-space formulation, this time due to cancellations arising from the rapidly os- 
cillating phase factor in the integrand in Eqn. (fTT?|) . We will occasionally distinguish the 
exponential suppression factors (which we can describe as 'dynamical' as they depend on 
the initial wave packet parameter Pq) from the more intrinsic pre-factor term (containing 
only the spreading time) which is due to the natural dispersion of the wave packet (which 
we can therefore describe as 'dispersive'.) 

For this case of an isolated quantum system, the specific result for the autocorrelation 
function in Eqn. ()2()|1 must also satisfy the general theorem in Eqn. (fT^ . To lowest non- 
trivial order in t ((9(t^)), the modulus of A{t) in this case for short times is 



|A(t)|2 = =exp 



^°(l + (t/2to^2^ 



{t/2tof 





" i-iir(''»+5?j+-- 

We note that this result arises from both the exponential ('dynamical') suppression, as well 
as the ('dispersive') prefactor term. The general result of Eqn. (fT^ requires the calculation 
of 



which give 



1 V 2 



When the right-hand-side of Eqn. (jl4j) is expanded to lowest order in t (again, 0{t'^)), using 
this result, it is found that 



fAHt\ f 1 fAHt\ 



2rn?a'^ 
~2t 




+ ■■■ (28) 



since to = rnha^ and the result in Eqn. is found to satisfy (in fact, to saturate, at this 
order) the inequahty in Eqn. ()14|) for short times, confirming the general result. We note 
that it is easy to show that the inequality is also satisfied to 0{t'^) where the left-hand-side 
is indeed larger than the right at this order. 



III. UNIFORM ACCELERATION 

The explicit form for Gaussian solutions, in both momentum-space and position-space, 
for the problem of a particle undergoing uniform acceleration (constant force F or linear 
potential given by V{x) = —Fx) are given in Ref. [s^. For the momentum-space form we 
have 



0(p, t) = $(p - _Pi)e-«P^V6m™ ^ g-a2({p-Ft)-po)V2 ^-i{p-Ft)xo/h ^i{{p-Ftf-p3)/emI^^Q^ 

(30) 



0(p,O) = J^e 



a 



-a^{p-poY /2 ^-ipxo/h 



with arbitrary initial position (xg) and momentum {po)- The corresponding position-space 
wavefunction is 



^iFt{xo-Ft^/6m)/h ^i{po+Ft){x-xo-pot/2m)/h 



1 



-Kah{l + it/to) 



X e 



-ix-{xo+pot/m+Ft^ /2m))^ /2{ah^)^{l+it/to) 



(31) 



with ■iIj{x,0) given in Eqn. (fTHj) . 

The calculation of A{t) can be done using either form to obtain 



A{t) 



exp 



{2tplt/mh - («Ft)2(l + (t/2to)2)) 



-iFt{xo-Ft'^/6m)/h 



(32) 



/l-it/2to ' L 4{l-tt/2to) 

and the same factors of 1 — it/2to as in Eqn. (fTn|l are obtained; this expression also reduces 
to that case in the free particle limit when F — 0, as it must. The modulus-squared is given 
by 



1 



^/TTJt/2to) 



exp 



-2a^(p^ + (Fto)ll + (t/2to) 



(t/2to) 



1 + (t/2to)' 



(33) 



and we note that this result can be obtained from Eqn. by the simple substitution 



Po 



pl+{Fto)\l + {t/2tor) 



(34) 



For this case of uniform acceleration, the wave packet spreading is identical (same Axf) as 
in the free-particle case, which can be understood by noting that the distance between two 
classical particles starting at the same initial location, undergoing the same force, but with 
slightly different initial velocities (or momenta, p\f^ — p^^^ = Apo) would be 

XAit) - xeit) = {xo + p[^h/m + Ff/2m) - (xq + pfh/m + Ft^ /2m) = (35) 

which increases linearly with time, in exactly the same way as for the free-particle solutions 
(when F = 0). The 'distance' between the peaks in ipo and ipt, however, eventually grows 
as t"^ so that the exponential ('dynamical') suppression in A{t) does not saturate, while the 
'dispersive' pre-factor is the same as for the free-particle case. 

We also note that the factors of po and F in |A(t)P appear in quadrature, and not in 
a combination such as po + Ft. One might naively expect that in cases where po and F 
have opposite signs, so that at a time given by tret = '^Po/ F when the classical particle (and 
central value of the quantum wave packet) has returned to Xq, the magnitudes of ilj{x,t) 
and ip{x,0) would be similar (only differing in possibly small spreading effects due to the 
[l + it/to) factors) and so would give rise to a relatively large value of |v4(t)|. While this does 
indeed occur for the magnitudes, at that time the classical momentum is of the opposite sign 
(p(0) = Po pitret) = —po), giving rise to rapidly oscillating phase factors in the quantum 
wave function, which still gives the expected exponential suppression. 



IV. SIMPLE HARMONIC OSCILLATOR WAVE PACKETS 

The free particle and case of uniform acceleration, corresponding to unbound motions, do 
not provide examples for direct comparison to the (quasi) periodic behavior seen in bound 
state systems. The harmonic oscillator admits Gaussian wave packet solutions which can be 
written (using propagator techniques) in closed form for arbitrary initial values of position 
and momentum {xq, po) and for which the evaluation of A{t) is therefore possible. In this 
case, we can also write any general wave packet as an expansion in eigenstates as in Eqn. (jSJ, 
using En = {n + l/2)huj, as 

oo 

^(x, t)=Y^ anWn(x)e-*("+i/2M (36) 

n=0 

from which it is clear that the observable probability density, |'?/'(x, t)p, is periodic with the 
classical period, = 2n/uj. 
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Using the propagator techniques outhned 
function 

^{x, 0) = 



1 



in js^, and the initial position-space wave 

(37) 



TT 



where (3 = ah, one can evaluate the time-development in closed form as 



1 



: exp 



where 



and 



L{t) = (3cos{ujt) 



S[x,t] 

/L(i)7?^^^" [2(3L{t)\ 
ih 



mujf3 



sin (cut) 



S[x,t] 



~xl cos(co't) + 2xxo — x^ 



imu}(5'^ sin(a;i:) 

h 



(38) 



(39) 



(40) 



2xoPo sin(co't) ^ 2ij3 pqx ijS pQsm.{ujt) 



The corresponding position-space probability density can be written as 



(x,t)p 



exp 



(x — Xq cos(u;t) — po sm{ut) / mu) 



(41) 



with 



{x)t = xo cos(ti;t) + 



Po sin(co't) 



and 



Axt 



\m\ 



(42) 
while the 



Thus, the expectation value moves in accordance with classical expectations 
width oscillates (from wide to narrow, or vice versa.) For the special case of the 'minimum 
uncertainty' wave packet where 



h 



rriLU 



the width of the packet is fixed as 



V2 



(43) 



(44) 



which is the same as the ground state oscillator energy eigenvalue state, but simply oscillates 
at the classical frequency. We note that the momentum-space wavefunctions can also be 
written using a propagator formalism Q and used to evaluate A{t) in a parallel fashion, 
obtaining the same result. 

The evaluation of |A(t)P for general values of /3, Xq, and Pq is straightforward enough, 
but the resulting expressions are somewhat cumbersome, so we will focus on several special 
cases as illustrative. 
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Case I: Minimum uncertainty wave packets, /3 = Pq. 
In this case the evaluation of A(t) gives 



A{t) = a/ cos(co't) + i sin(u;t) exp 



X, 



_|_ (^oPo 



2ff 



[(1 — cos(c<jt)) — i sin(co't)] 



where great simphfications have been made by noting that 

cos (cut) + i sin (cut) 



(45) 



(46) 



cos (cut) — i sin (cut) 

We note that a very similar expression arises in analyses of the macroscopic wavefunction 
for Bose-Einstein condensates {sS] and the colla pse and revival of the matter wave field for 
such systems has been observed experimentally j3J]. 

Once again, the two important parameters appear together in quadrature, as in the 
uniform acceleration case. This gives 



\Ait)\' 



exp 



X, 







) 



[1 — cos(ct;t)] 



(47) 



which clearly exhibits the expected periodicity. All of the suppression can be attributed 
to the 'dynamical' factors (those in the exponential, containing xq and po) as there is no 
'dispersive' pre-factor component for this constant width packet. 

For this case, the minimum degree of overlap at any point during a single classical period 

is 



\A{TJ2)\' 



exp 



-2 ( ^ + 



(48) 



so there is no time at which the wave packet is ever truly orthogonal to its initial state. 



Case II: Arbitrary /3, but Xo,Po = 0. For this case, the wave packet does not oscillate, but 
only 'pulsates', and the time- dependent wave function simplifies to 



1 

ip{x,t) = ^ exp 

i/vry {(3 cos{ujt) + {ih/ muojS) sin(ci;t)) 



— x^[cos(c<jt) + (imuP'^/h) sin{ujt)] 



2/?2[cos(cut) + {ih/muj(3) sin(cut)] 

(49) 

It is convenient to define the parameters 

r = ^ = 9 so that - = ^ 50 

in terms of which the resulting autocorrelation function in this case has the very simple form 



^' 2cos(cut) -i(r + l/r)sin(cut) ^^^^ 
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or 



\A(t)\^ = ^ (52) 
v/cos2(cut) + (r + l/r)2 sm\ujt)/4 

all of which can be attributed to a 'dispersive' (but in this case periodic) pre-factor. 

We first note that in this case A{t) is invariant under the transformation r ^ 1/r, in other 
words, the time-dependence is the same for both initially wide (/3 > Po) or narrow (/5 < Pq) 
packets. It is clear that the larger the deviation from the 'minimum uncertainty' wavepacket, 
the faster the wavepacket 'pulsates' away from its initial shape. It is also noteworthy that in 
this case |A(Tc;/2)| = 1 so that the wave packet returns to its initial form (up to a constant 
complex phase) twice each classical period. This can be understood from the expansion of 
this wave form in terms of energy eigenstates. In this case, where the parameters Xq,Po both 
vanish, one is expanding an even-parity function in Eqn. (jSE)), so that only the even (a2n) 
terms are nonvanishing and the n-dependent exponential factors in Eqn. ^ oscillate twice 
as rapidly as in the general case. 

Finally, for the very special case where (3 = /Sq [r = 1) as well, we recover the ground 
state energy eigenstate of the oscillator, with its trivial stationary-state time-dependence 
{ipo{x,t) = uo{x) exp{—iEQt/h)) and Eqn. indeed reduces to 



A{t) ^ J- , , \. . , , = = e+*-*/2 (53) 

^ ^ Y 2 cos{tut) - 2i sm{ujt) ^ ' 

as expected. 

Wave packet solutions (Gaussian or not) of the harmonic oscillator can be shown (using 

n 

the expansion in Eqn. (jSEI), for example) to satisfy [2^ 

^(a;,t + mr,0 = (-l)'"^(x,t) (54) 

with a similar result for the momentum-space version as well. Because of the specially 
symmetric nature of the potential, we also have 

V'(-x, t + T,;/2) = (-i) ^(x, t) and (/.(-p, t + T,,/2) = {-%) (f){p, t) (55) 

so that half a period later, the wave-packet is reproduced, but at the opposite 'corner' of 
phase space, namely with x ^ —x and p ^ —p: we note that two applications of Eqn. ()55|1 
reproduce Eqn. (jJlj). One can also show these connections using the propagator techniques 



in Ref. 
in Ref. 



3C|, provided one properly identifies the complex pre-factors as described in detail 
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This type of behavior can be diagnosed using a variation on the standard autocorrelation 
function, namely 

r" + 00 



iIj*{—x, t) ip{x, 0) dx 



(f)*{-p,t)<j){p, 0)dp 



(56) 



which measures the overlap of the initial state with the 'out-of-phase' version of itself at 
later times. Given the simple connections in Eqn. (fS^ . we can immediately write, for Case 
1 considered above, 



exp 



+ COs(co't) 



(57) 



which is exponentially suppressed at integral multiples of T^, but unity at t = {2k + l)Tc//2. 
This type of anti- correlation function finds use in the study of wave packet revivals ^ \ where 
quantum wave packets may reform near t = T^ev, as in Eqn. (jl2|) . but out of phase with the 
original packet. 



V. 'INVERTED' OSCILLATOR WAVE PACKETS FOR UNSTABLE EQUILIB- 
RIUM 



As described in detail in Ref. j30|], the case of the 'inverted' oscillator, defined by 

V{x) = —niuj'^x^ (58) 

can be studied using wave packet results for the standard harmonic oscillator by making the 
substitutions 

u ^ iu , sm^Lut) ismh.{o)t) , and cos(u;t) cosh((Dt) . (59) 

The general wave packet solution in Eqn. ()38p. for example, can be carried over in this way 
to obtain the 'runaway' wavepacket, with probability density given by 

(x — Xo cosh(a;i:) — Pq smh{ut) / moj)^ 



V^\B{t)\ 



(60) 



with 



where 



, V , /- X Posinhfcut) , ^ \B(t)\ . 

{x)t = xocosh(ujt) + — ^ and Ax* = (61) 



\B{t)\ = \/(3^ cosh^ (Cot) + {h/mu(3) ^ sinh^ (ut) . (62) 
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As above, the expression for A{t) for the general case is cumbersome, so we only examine 
it for one specific case as an example, namely the case where (3 = (3^ = \/ h/muj. This 
situation no longer corresponds to a constant width wave packet, since 



Ax. 



cosh {ijjt) + sinh {ujt) 



(63) 



increases exponentially, as the individual momentum components comprising the wave 
packet quickly diverge in p-space. For the case of xq = 0, we have the general expression 



A{t) 



1 



y/cosh^cut) 



exp 



2mujh 



cosh(co't) — 1 + z sinh(ti;t)(2 cosh(ci;t) — 1) 



(64) 



cosh (cut) (cosh (cut) — 2 sinh (cut)) 

In the limit when t » 1/uJ, the hyperbolic functions both approach exp (cut) /2 and we have 
the limiting case 

1 



A{t » 1/Cd) 



exp 



2muh 



(65) 



The exponential ('dynamical') suppression once again is seen to 'saturate', as in the free- 
particle case, and for the same reason, namely that both x{t) — xq and Axt have the same 
large t (here exponential) behavior. The resulting modulus is given by 



2e~^^ exp 



mujh 



(66) 



which still becomes exponentially small, but now due to the ('dispersive') prefactor. If one 
also has xq ^ 0, the expression above includes an additional factor of exp(— Xq//?q) (similar 
to that in Eqn. (jTfjl . with no cross-term involving xq times po- 



VI. CONCLUSIONS AND DISCUSSION 



We have evaluated the autocorrelation function for time-dependent Gaussian wave packet 
solutions for four quite different quantum mechanical systems, examining the behavior of 
A{t) in terms of both classical analogs and the quantum mechanical time-evolution of the 
wave function magnitude and phases in both position- and momentum-space. We have 
focused attention on the different contributions to the suppression in A{t) arising from 'dis- 
persive' (prefactor) and more 'dynamical' (typically exponential) effects (which depend on 
the initial conditions) which may or may not saturate to a small but finite value, depending 
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on the relationship between the classical dynamical behavior {x(t) or {x)t) and the quan- 
tum mechanical spreading (Axt). We have extended earlier results of Baltz in terms 
of a now standard analysis tool, namely the autocorrelation function, directly comparing 
the behavior of A{t) for several distinct classes of classical behavior. Using the free-particle 
solution, we have also been able to exhibit a useful test case for quite general theorems on 
the time-development of isolated quantum systems, while providing other non-trivial exam- 
ples of closed-form results for A{t) for Gaussian solutions. Additional examples which are 
simple extensions of the results presented here are multi-dimensional free-particle, uniformly 
accelerated particle, or harmonic oscillator solutions, where the autocorrelation function fac- 
torizes as A{t) = Ax{t) ■ Ay{t) , or the related problem of a Gaussian wave packet in a uniform 
magnetic field which admits Gaussian solutions corresponding to classical circular orbits. 
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